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We propose a new method to constrain elastic scattering between dark matter (DM) and standard 
model particles in the early Universe. Direct or indirect thermal coupling of non-relativistic DM 
with photons leads to a heat sink for the latter. This results in spectral distortions of the cosmic 
microwave background (CMB), the amplitnde of which can be as large as a few times the DM-to- 
photon number ratio. We compute CMB spectral distortions due to DM-proton, DM-electron 
and DM-photon scattering for generic energy-dependent cross sections and DM mass > 1 
keV. Using FIRAS measurements we set constraints on the cross sections for 0.1 MeV. 

In particular, for energy-independent scattering we obtain rroM-proton < 10“^^cm^ (keV, 
CTDM-eiectron < 10“^^cm^ {ksN and CTDM-photon < 10“®®cm^ (mx/keV). An experiment with 
the characteristics of PIXIE would extend the regime of sensitivity up to masses ~ 1 GeV. 


Introduction - In the standard cosmological model, 
dark matter (DM) only interacts gravitationally, and a 
single number (its abundance) suffices to accurately de¬ 
scribe its effect on the cosmic microwave background 
(CMB) anisotropies and large-scale structure (LSS). Yet, 
a variety of DM candidates are predicted to interact 
weakly with standard model particles mm- Several di¬ 
rect detection experiments HH, and various astronomi¬ 
cal and cosmological probes mm have or will set limits 
on interacting DM models, but large portions of param¬ 
eter space are still unconstrained. 

In this Letter, we propose a new probe of DM interac¬ 
tions with standard model particles, relying on spectral 
distortions (hereafter, SDs) of the CMB blackbody spec¬ 
trum. It is well known that energy injection at redshifts 
0 < 2 X 10® can distort the CMB spectrum [51110]. SDs 
take the form of a chemical-potential (“/r-type”) distor¬ 
tion if energy injection occurs early enough (x > 5 x 10^) 
that Compton scattering efficiently redistributes photons 
in frequency. The SD shape smoothly transitions to a 
“y-type” distortion at lower redshifts when photon fre¬ 
quency redistribution becomes less efficient. SDs have 
been used, for example, to constrain decaying particles 
[n EH or the primordial small-scale power spectrum 
mM- Conversely, any energy extraction from the pho¬ 
tons can also induce SDs. Such a process is at work in 
the standard cosmological scenario EMI]: CMB pho¬ 
tons heat up the non-relativistic plasma through Comp¬ 
ton scattering followed by rapid Coulomb interactions. 
Instead of cooling down adiabatically with a temperature 
Th oc 1/a^, where a is the scale factor, electrons, protons 
and helium nuclei (“baryons” for short) are maintained 
in thermal equilibrium with the photons down to redshift 
z « 200, with Tf, « r.y oc 1/a jlH]- This process removes 
heat from the CMB and leads to negative /r-type and 
y-type distortions of a few times the baryon-to-photon 
number ratio, of order a few parts in a billion. 

The effect we describe in this work is an extension of 
the aforementioned process. If a non-relativistic DM par¬ 


ticle is thermally coupled to the electron-nucleon plasma 
through frequent collisions, the energy extracted from 
the photons by Compton scattering is then redistributed 
among a larger number of particles. Maintaining the 
baryons and the DM in thermal equilibrium with the 
radiation therefore requires an increased rate of energy 
extraction from the photons. Direct scattering of DM 
particles with photons have the same outcome. In both 
cases SDs are generated, of order up to a few times the 
DM-to-photon number ratio. Given the known DM mass 
density, this number is inversely proportional to the DM 
mass. 

To be more precise, this effect only takes place as long 
as the heating of the DM is more efficient than the com¬ 
peting adiabatic cooling due to cosmological expansion. 
Hence the final SD is at most a few times the DM-to- 
photon number ratio. It is smaller the shorter the epoch 
of tight thermal coupling of DM with the photon-baryon 
plasma is. For a given sensitivity to SDs, there is there¬ 
fore a maximum mass that can be probed, such that the 
SD reaches the instrument sensitivity when the DM is 
tightly coupled at all relevant times. This maximum mass 
is ^ 0.1 MeV for FIRAS [12] and will be ~ 1 GeV for a 
future SD experiment with the sensitivity of PIXIE [20| . 
Below that maximal mass, SD measurements can set an 
upper bound to the duration of tight coupling, hence to 
the scattering cross section. 

Since the maximal SD is inversely proportional to the 
DM mass, the effect we introduce allows to test light- 
DM models, while most direct-detection experiments are 
insensitive to sub-GeV masses. In this first study we 
limit ourselves to DM masses greater than a keV, such 
that the DM is already non-relativistic at the beginning 
of the distortion epoch, z « 2 x 10®. 

DM scattering with protons or electrons - We 

assume that the DM particle y can elastically scatter off 
baryons with a momentum-transfer cross section of the 
form a{v) = anv'^, where n is an arbitrary integer and v 
is the magnitude of the baryon-DM relative velocity (in 


2 


units of the speed of light). Some of the best-motivated 
DM models indeed have interactions of this form [2T] . For 
instance, electric-dipole or magnetic-dipole DM coupling 
to the standard model through heavy charged messengers 
have cross sections a (x and v'^, respectively. In this 
Letter we specifically consider scattering with either pro¬ 
tons or free electrons, and neglect scattering with helium 
nuclei. 

Baryons maintain a single temperature at all times 
due to extremely frequent Coulomb scatterings [18] . The 
baryon temperature is in turn closely coupled to the 
CMB temperature through Compton scattering of 
photons by free electrons. As long as the rate of DM- 
baryon scattering is much larger than the Hubble expan¬ 
sion rate, DM particles have a Maxwellian velocity distri¬ 
bution with temperature ^ Th. Once the collision rate 
falls below the expansion rate, then {i) the DM momenta 
start redshifting freely as oc (l-|-z), and {ii) if this de¬ 
coupling is not instantaneous and is velocity-dependent, 
the DM velocity distribution is no longer nec¬ 

essarily described by a Maxwell-Boltzmann law. One 
should then in principle compute by solving the Boltz¬ 
mann equation. 

In the limit that the DM distribution is thermal and 
that random velocities dominate over bulk flows (valid 
for z > 10'^), Ref. [52] have showed that the evolution of 
the DM temperature is governed by 


with 


t^ = -2HT^ + T^,{T,-T^), ( 1 ) 

_ 2cnNi,a„mi,m^ / Tj 
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where c„ is a constant depending on n given in [22] 
and Nb = is the number density of scatter¬ 

ing baryons. Assuming radiation domination (valid at 
z > 3000), the Hubble expansion rate H{a) scales as 
H = Setting Tb = the ratio 

r^b/H takes the form 


Fxb _ ( m^/mb + T^/Tb \ " 

H \ a ) \ m^/rub + l ) ’ ^ ’ 


where we have defined the characteristic scale factor a^^b 
such that: 

The scale factor a^b marks the transition between ther¬ 
mal coupling (Fjj-b > H) and decoupling (F^;, < H). In 
this work, we will only consider slopes n > —3 for which 
the DM starts tightly coupled to the baryons at early 
times, and thermally decouples at a > a^b- In these cases 
DM-baryon scattering can have negligible effects on CMB 
anisotropies and LSS, yet can still manifest itself at high 
redshifts and induce SDs. 

We define the dimensionless parameter 



characterizing the efficiency of the DM-baryon thermal 
coupling. In the limit that Fj^-f, ^ iL, DM and baryon 
temperatures are tightly coupled, = Ty, which 

implies « —HT^. Inserting this value back into the 
left-hand-side of Eq. § we obtain r^b{oi ^ o-xb) — 1- On 
the other hand, r^bi^ ^ axb) 0 since Fy^,;, ^ iL in that 
regime. The transition from tight coupling to decoupling 
is sharper as n is increased, as can be seen from Eq. 

Eor n = —2, T^b/H ~ y/a^pj~a and thermal decoupling 
is very gradual, taking a few decades in scale factor. This 
implies that the description of the DM velocity distribu¬ 
tion with a Maxwell-Boltzmann law is very inaccurate 
in that case (unless the DM is self-interacting, which we 
do not consider here). We defer the full solution of the 
collisional Boltzmann equation to future work, and will 
focus on n > — 1 in what follows. 

We confirm all these features by solving Eq. Q nu¬ 
merically and computing the resulting r^f,(a). We find 
that r^b ~ 1/2 at a « {2/3)a^b for a broad range of 
mass ratios and slopes n > — 1. Once the DM thermally 
decouples from the baryons, its velocity distribution is 
no longer Maxwell-Boltzmann and Eq. Q is no longer 
valid. As a result the parameter r^b(a} we computed is 
inaccurate for a > a^-b- We therefore adopt the follow¬ 
ing simple step-function approximation for r^bia), which 
captures its essential features: 


r^bio-) = 1 for a < (2/3)a^t, and 0 otherwise. (6) 


DM-photon scattering - We also consider direct 
scattering of DM particles with photons, with an energy- 
dependent momentum-transfer cross-section a'(E^) = 
ap{Ey/E qY. We set the normalization at Eq = 1 
keV, close to the characteristic CMB photon energy at 
z ~ 10®. A Thomson-like scattering {p = 0) would arise, 
e.g., for a millicharged DM. A quadratic dependence on 
energy (p = 2) could arise from a DM particle with an 
electric or magnetic dipole moment [53]. A quartic de¬ 
pendence {p = 4) occurs e.g. for Rayleigh DM [53]. As¬ 
suming the DM is non-relativistic {m^ 3> I/), the result¬ 
ing heating rate is easily obtained by generalizing the cal¬ 
culation of the Compton heating rate due to scattering 
off free electrons [55] : 


^xlx7 “ 


with 
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(7) 

( 8 ) 


where dp is a numerical constant, with value (0.28, 1, 4.8, 
28.2, 1558) for p = (—1,0,1,2,4) respectively. 

In this case, the ratio of interaction to expansion rates 
takes the form 


H 

where (a^.y)^'*'^ 




8 dpapp° 


{T^/EoY. (10) 
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We only consider p > — 2 so that DM-photon interac¬ 
tions are efficient at early times and irrelevant at a > 
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cLxt Here again we define the dimensionless parameter 
~ '^x)l compute it numerically 
(assuming DM scatters off photons only). Its behavior is 
very similar to the one obtained for rxb^ with the corre¬ 
spondence n e-)- 2p -I- 1, and we also approximate it by a 
step function with a transition at a = {2/3)ax'y 

Spectral distortions - The baryon temperature 
evolves according to (i) adiabatic cooling due to cosmo¬ 
logical expansion, (ii) Compton heating by CMB photons 
and (in) energy exchange with the DM: 

n = -2HT, + TcoAT, - n) + -^Txb(Tx - n). (11) 


Here = Px/iTix is the number density of DM particles 
and Nl°^ is the total number density of “baryons” (nuclei 
and free electrons) maintained in equilibrium with one 
another. As long as the plasma is fully ionized (valid for 
0 > 6000), Nl°^ = 2Nn+iNne « p6/mH(2-|rHe) where 
A^h and A^ue are the abundances of hydrogen and helium 
and lije is the helium fraction by mass. The last term in 
Eq. (11) is easily obtained from Eq. Q by requiring that 
baryon-DM collisions conserve the total thermal energy. 

The rate of extraction of energy from the photons by 
Compton scattering is [T5] 


Px 


dt 


Py 


Com 


= -NI°^TcoM - Ty). (12) 


Since for all redshifts z > 200 Compton scattering main¬ 
tains Tb « Ty to very high accuracy, the rate of change 
of the baryon temperature is just Tb = —HTb- Replacing 
the left-hand-side of Eq. 0 by this value we obtain the 
net Compton heating rate, and substitute it in Eq. (12) 
to arrive at 


= -l{Nl°^+rxbNx)HT,, (13) 
ai \ pj / z 

where the parameter was defined in Eq. (§. The first 
term was derived in Ref. m and translates the extrac¬ 
tion of energy from photons due to Compton heating of 
the baryons. The second term arises from the enhanced 
heat capacity of the baryon-DM fluid due to DM-baryon 
scattering. 

The effect of direct DM-photon scattering is similar: 
in this case the cooling rate of the photons is given by 

P^Jt = -ifxyNxHT^- (14) 


[B M- Compton scattering by free electrons in hot 
clusters and the reionized intergalactic medium leads to 
a y-distortion of ^ 2 x 10“® [57]. This is below the sen¬ 
sitivity of FIRAS [ini, but orders of magnitude larger 
than that of PIXIE . We shall therefore not consider 
pure ^-distortions in this work, and cut the integration 
at Zmin = 10"^ (flmax = 10“'^). This allows us to con¬ 
sider a radiation-dominated and fully ionized universe, 
and neglect bulk flows relative to thermal velocities (55) . 

Our final estimate for the amplitude of SDs due to 
DM-baryon collisions is therefore 


A = 


Ap, 
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where A^ « pj/(2.1T^) is the number density of 
CMB photons, and we have taken the constant baryon- 
to-photon and DM-to-photon number ratios out of 
the integrals. The cutoff a* is defined as a* = 
max[ayi, min(aniax; (2/3)0^^;,)] (we smooth the transitions 
for better visual results). The same expression applies to 
DM-photon scattering with the substitution axb dxy 
The integrals only depend logarithmically on the bound¬ 
aries. Our various approximations (i.e., taking a sharp 
boundary at = z~^ instead of using a SD visibility 
function m, assuming a step function for rxb(a) and 
choosing Umax = 10“^) should therefore not significantly 
affect our results. They have the advantage of giving 
simple analytic expressions. 

Results - All our results are computed with the cur¬ 
rent best-fit values for cosmological parameters [HI [28] . 
We show in Fig. the photon distortion ^p^jp^ for 
velocity-independent DM-proton scattering as a function 
of (To, for several values of the DM mass. For cross sec¬ 
tions small enough that (2/3)0;^;, < (a* = a^j), DM 

scattering has no effect and the distortion plateaus at 
Ag « —3x10“®, due exclusively to the cooling of baryons 
[IKHTT] . For cross sections large enough that (2/3)a^6 > 
Omax, the DM is tightly coupled to baryons at all rele¬ 
vant times, and the distortion saturates (a* = Omax re¬ 
gardless of the cross section). In this limit, the distortion 
is enhanced by a factor A^/A^°* « 3(GeV/m^), so that 
A[(t„ —)• oo] « Ag — 9 X 10“®(MeV/m^). Between these 
two regimes, the distortion scales logarithmically with 
Gxb and hence (t„. 

For a given sensitivity Amax, a maximum mass 
can be probed, such that |A[(t„ —\ = A^ax: 


To obtain the final relative amplitude of spectral dis¬ 
tortions, we have to integrate Eq. (13) or (14) over time. 
The high-redshift boundary is at « 2 x 10®: energy 
injection (or extraction) at z > simply leads to a 
change of the photon temperature and does not distort 
the blackbody spectrum [10]. The low-redshift end is 
in principle Zmin « 200, the epoch of thermal decou¬ 
pling of baryons from photons. In practice, energy in¬ 
jection at z < 10"^ leads mostly to a y-type distortion 


m 


max 

X 


1 MeV X 


9 X 10“® 

A„,ax - 3 X 10-9 ■ 


(16) 


Higher masses are completely unconstrained as they lead 
to a too small number density of DM particles. FIRAS 
has measured the CMB blackbody spectrum to an ac¬ 
curacy Ainax « 5 X 10“® [H], implying « 0.18 

MeV. The proposed experiment PIXIE [50] will reach a 
sensitivity Amax ~ 10“®. This will allow to constrain 
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FIG. 1. Absolute value of the photon distortion for 

DM collisions with protons, for a velocity-independent cross 
section (tq. The solid curves are labelled by the DM particle 
mass. The upper dashed curve indicates the approximate 
constraint from FIRAS Ap-y/p-y < 5 x 10“® |19) . The lower 
dotted curve indicates the approximate forecasted sensitivity 
of PIXIE Ap-^/p^ ~ 10“® [ig. 


the DM-baryon and DM-photon cross sections up to DM 


masses m' 


1.3 GeV. 


For masses rriy^ < m™®’’', 


a measurement of the CMB 
blackbody spectrum to a precision A^ax would imply an 
upper limit on the cross sections cr^^ o r por DM- 

baryon collisions we obtain, using Eqs. (15) and (|^, 




mt 


mb 


my, 




.(17) 


For DM-proton collisions, the numerical constants Cn are 
(1.4 X 10-30,1.1 X 10-27,8.2 X 10-20,5.5 x 10-22) cm2 
for n = (—1,0,1,2) respectively. For DM-electron col¬ 
lisions, the corresponding values are (1.4 x 10-30,2.6 x 
10-20,4.5 X 10-23,7.0 X 10-27) cm2, constraint on 
the DM-photon cross section is obtained similarly from 
Eqs. (15) and (10): 


P ^ ^MeV 


(P+2)mx/"i; 


(18) 


with Dp = (6.3, 5.6, 3.7, 2.0, 0.4) x 10-37 cm2 £qj. ^ _ 
(—1,0,1, 2,4), respectively. 

Equations and ([T^ are the main results of 


this Letter. Given a sensitivity Amaxj they allow to ob¬ 
tain upper limits on DM-baryon and DM-photon cross 
sections with power-law dependence on the baryon-DM 
relative velocity or photon energy (with n,p > —1), up 
to a maximal DM mass 

We plot in Fig.j^the current constraints on the energy- 
independent cross sections ctq^, ctq'^, CTq^' as a function 
of the DM mass given the FIRAS measurements. We 
also show the forecasted constraints for the sensitivity of 
PIXIE. 

Comparison with previous bounds ~ Most direct 
detection experiments only constrain DM-nucleon cross 


FIG. 2. Current upper bounds from FIRAS (solid) and 
forecasted detection thresholds from PIXIE (dotted) on the 
energy-independent DM-proton (purple), DM-electron (blue) 
and DM-photon (red) cross sections no, as a function of the 
DM mass. Masses > 0.18 MeV are unconstrained by FI¬ 
RAS as the distortion can never reach Ap-y/p^, = 5 x IQ-®, 
even for infinititely large cross section. PIXIE will extend 
the domain of constrainable masses by four orders of mag¬ 
nitude, up to nix ~ 1.3 GeV. For comparison, we also show 
the constraints on DM-electron scattering from XENONIO 
data [6] and the limits on DM-photon scattering from Milky 
Way satellite counts . No other probe currently constrains 
DM-proton scattering in the range of masses shown. 


sections for masses > few GeV, required to produce 
sufficient nuclear recoil. Ref. |22] derive constraints on 
the ratio for DM-proton collisions in the limit 

mn, using GMB anisotropy and LSS data. SDs 
therefore provide a probe of DM-nuclei scattering in a 
mass range complementary to the one currently con¬ 
strained. In particular, our limits on DM-proton scat¬ 
tering from FIRAS measurements are the only existing 
bounds for rriy^ ($0.1 MeV. 

Ref. [S] have set the first constraints on the scattering 
of sub-GeV DM with electrons, which could lead to ion¬ 
ization events in the target material |30j . For a velocity- 
independent cross section, they find co $ 3 x 10-3® cm2 
for rriy^ = 100 MeV, significantly better than what we 
forecast at the same mass for a PIXIE-type experiment, 
CTo $ 10-23 cm2, bound of Ref. [5], however, worsens 
rapidly for DM masses below a few MeV. Here again, FI¬ 
RAS limits give the only existing bounds on DM-electron 
cross sections for rriy^ $ 0.1 MeV. 

Ref. |53] give a constraint on the DM-photon energy- 
independent cross section using counts of Milky Way 
satellites, translating to < 3.7 x 10-33(TO^/MeV) cm2. 
The constraint we set with FIRAS for my^ <C 0.1 MeV is 
tighter by a factor of ^ 5, and PIXIE will allow to ex¬ 
tend it up to my^ « 1 GeV. We also constrain the p = 2 
cross section (72 $ 2 x 10-37(TOjj./MeV), tighter by six or¬ 
ders of magnitude than the limit of Ref. [3T] using GMB 
anisotropies. 

Conclusions - We have set forth a new avenue to 
probe DM interactions with standard model particles. 
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using CMB SDs. We have studied the effect of DM scat¬ 
tering with either protons, electrons or photons, for a 
power-law velocity and energy dependence of the cross 
section. We have shown that the FIRAS measurements 
can already set constraints on the cross sections for DM 
masses < 0.1 MeV. Above this mass, the number 
of DM particles is too small to affect the effective heat 
capacity of the plasma at a sufficient level. The high 
sensitivity of PIXIE will allow to constrain DM particles 
with masses up to ^ 1 GeV. 

Specific models for the DM particle would predict the 
shape and relative strengths of interactions with differ¬ 
ent species. The overall SD can be simply obtained by 
linearly adding the contributions of each scattering pro¬ 
cess. While we have only considered the characteristic 
amplitude of the distortion in this work, its detailed form 
can easily be obtained by convolving the energy extrac¬ 
tion rate with a SD Green’s function [26| . In particular. 


when DM-baryon or DM-photon decoupling occurs after 
z ^ 5 X 10^, SDs may allow determining specific DM pa¬ 
rameters through the residual (non-p, and non-j/) distor¬ 
tion [13 [MISS]- Delving in such details would, however, 
require a better treatment of the DM velocity distribu¬ 
tion, which should be computed by solving the Bolzmann 
equation rather than assuming it is Maxwellian. We have 
also only considered DM masses > 1 keV, such that 
the DM particle is non-relativistic at z < 2 x 10®. It 
would be interesting to extend this work to lower masses 
and relativistic DM particles. We defer the study of these 
questions to future work. 
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